Correlations in an expanding gas of hard-core bosons 
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We consider a longitudinal expansion of a one-dimensional gas of hard-core bosons suddenly 
released from a trap. We show that the broken translational invariance in the initial state of the 
system is encoded in correlations between the bosonic occupation numbers in the momentum space. 
The correlations are protected by the integrability and exhibit no relaxation during the expansion. 
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Rapid progress in the ability to manipulate ultracold 
atomic gases stimulated a revival of interest in funda- 
mental properties of interacting Bose systems. Unlike 
conventional condensed matter systems, cold gases of- 
fer a unique possibility to monitor the out-of-equilibrium 
dynamics of interacting systems unhindered by coupling 
to the environment and associated with it decoherence, 
sec [1] for a recent review. Among various realizations 
of interacting Bose systems, the one-dimensional (ID) 
ones @, H, 0| occupy a special place: interactions in ID 
have a much stronger effect than in higher dimensions 
while often allowing for a complete theoretical treatment. 

This paper is partially motivated by the recent experi- 
mental study Q of the evolution of ID strongly interact- 
ing Bose liquid from a carefully prepared nonequilibrium 
initial state. The experiment [3] showed that the momen- 
tum distribution function does not exhibit a noticeable 
relaxation towards equilibrium. Such extremely slow re- 
laxation is consistent with the behavior expected for an 
almost integrable system. Indeed, in a ID system only 
three-particle collisions may lead to a momentum relax- 
ation; such processes are absent in integrable models 

On the theoretical side, experiments such as 0,01 high- 
light the relevance of the well-known in statistical me- 
chanics quantum quench problem: how to describe the 
evolution of a system from an arbitrary initial state (see, 
e.g., and references therein). At present, there are 
very few exact results on such strongly nonequilibrium 
dynamics of interacting quantum systems. Various prob- 
lems of this type arise naturally in the description of the 
experiments on trapped cold atomic gases. Indeed, by 
far the most popular technique today is to observe the 
expansion of a gas after a sudden release of the trap [l| . 
Although such experiments are obviously destructive, the 
time-of-flight imaging [l[ allows one to study the real- 
time evolution of the bosonic occupation numbers in the 
momentum space. Importantly, not only the average oc- 
cupation numbers (momentum distribution) but also the 
corresponding higher-order statistical moments (fluctua- 
tions) are accessible experimentally [7|, |8[). Unlike the 
momentum distribution, the fluctuations are sensitive to 
the relaxation in the system. 

We consider a simple yet realistic [2] example: expan- 



sion of a ID gas of bosons with infinitely strong con- 
tact repulsion (hard-core bosons) suddenly released from 
a trap. We show that shortly after the trap release, 
bosonic occupation numbers reach their steady-state val- 
ues. We derive an operator identity, see Eq. (p?0)) be- 
low, that relates the bosonic occupation numbers in the 
steady state to the integrals of motion. The identity al- 
lows one to study all statistical moments of the bosonic 
occupation numbers, and we evaluate the second moment 
in a closed form. Correlations between the occupation 
numbers at different momenta reflect directly the broken 
translational invariance in the initial (trapped) state of 
the system. 

To be specific, we assume that initially (at t < 0) the 
system is in a thermal equilibrium state of the Hamilto- 



H = Hn 



Vt 



trap ■ 



(i) 



Here Ho describes ID hard-core Bose gas without con- 
finement (see Eq. (j4|) below) and 



V t 



trap 



dx V(x)p(x) 



(2) 



describes the effect of the trap, with p(x) being the local 
density operator. The trap potential ((2]) breaks trans- 
lational invariance, hence it does not commute with the 
Hamiltonian Hq that governs the dynamics after the trap 
release at t = 0. 

The simplest description of the hard-core bosons is 
based on the Jordan- Wigner transformation 



exp 



dyp{y) 



(3) 



where the operators ip and <p correspond to fermions and 
bosons, respectively: 

{i/j(x),tp i (y)} = [(f(x), ^(y)] = S(x - y). 

The transformation ([3l) maps the hard-core bosons onto 
free spinless fermions 0], 



2m dx 2 



tp(x) 



(4) 
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At the same time, the local density operator [and, there- 
fore, Eq. ([2J], retains its form, 



p(x) = ip (x)tfi(x) = ip'(x)<p(x). 

Since the fermionic occupation numbers in the momen- 
tum space 

n k =4i>k, A = (^r 1/2 fdxe~ lh *iP(x), (5) 



commute with Hq, the expectation values of n k are inde- 
pendent of time, 



(n k ) f = (n k ) = const. 



(0) 



Hereinafter 



(0) t = (e ,Hti Oe- ,fl » f ) , 



where (. . .)o = (•••)t^o denotes the thermal averaging 
with the initial Hamiltonian H, see Eq. ((T|). 

The crucial for the following observation is that even 
though the evolution of the system at t > is governed by 
the translationally invariant Hamiltonian Hq, the initial 
Hamiltonian H does not have this symmetry. Therefore, 
not only the diagonal in k fermionic bilinears (such as n k ) 
have finite expectation values, but also the off-diagonal 
ones, e.g., (ifi k ifi k ,)o 0- Since the expectation value 



: p 2 /2m, (7) 



oscillates with t, quantities such as | (V'fc^/v)* | remain 
constant and carry with them the memory of the broken- 
symmetry initial state of the system. Because the off- 
diagonal correlation functions ([7]) are finite, the fermionic 
occupation numbers fluctuate. Indeed, with the help of 
the Wick theorem one finds for Sn k = n k — (n k ) 



(5n k Sn k ,) t = -j^^ol 



const. 



(8) 



It should be emphasized that the very survival of quan- 
tities such as Eqs. {7J and ||5J) that preserve the infor- 
mation about the initial conditions does not rely on the 
particularly simple form that the Hamiltonian H has in 
our case. Rather, it is a direct consequence of the in- 
tegrability. Indeed, in a generic (nonintegrable) system 
correlation function {7} would decay with t. This decay 
(relaxation) "washes out" the memory about the symme- 
try of the initial state, thereby restoring the translational 
invariance. After the relaxation is complete, the density 
matrix commutes with the total momentum. Averaging 
with any density matrix that has this symmetry would 
give zero for the correlation functions such as Eq. ([7]). 

A very similar consideration can be applied to any 
quantum quench problem in which the symmetry of the 
Hamiltonian that governs the system's dynamics differs 
from that of the initial state. [For example, suddenly 



turned off interactions in the Luttinger model [10( corre- 
spond to the initial state with broken global U(l) symme- 
try] . The information about the symmetry of the initial 
state is encoded in the off-diagonal correlation functions 
[cf. Eq. 0]; relaxation manifests itself in the decay of 
these off-diagonal correlations with time. 

We now consider a specific but rather realistic situation 
when the trap potential Eq. ([2]) is harmonic, 



2 2ml 



(9) 



At zero temperature the correlation function ((7|) can be 
written as 



N-l 



(10) 



where N 3> 1 is the number of particles in the system 
and <fi n (k) is the stationary eigenfunction in the momen- 
tum representation that corresponds to n-th energy level 
of a harmonic oscillator. The expectation values of the 
fermionic occupation numbers ([6]) are obtained by setting 
k = k' in Eq. (Tl0|) . For N ^> 1, this yields a "semicircle" 
dependence 



(n k ) = -Jl-k*/k* F , 

7T » 



(11) 



where kp is the Fermi momentum and R is the classical 
radius of iV-particlc fermionic cloud confined in a har- 
monic trap, 

k F l = R/l = V2N. 

In writing Eq. (fTTj) we neglected the oscillating with k 
contribution that has a relative magnitude of the order 
of l/R 1. This contribution is the momentum-space 
counterpart of the Friedel oscillations in (p(x)) , as it is 
obvious from the operator identity 



n k 



AHt 



[l 2 p(kf 



7T 

2lo 



(12) 



The period of the Friedel oscillations in x-space is the 
Fermi wavelength 2-K/kp- Eq. (fT2|) then implies that 
the corresponding oscillations in fc-space have a period 
2tt/M 2 = 2k/R- 

Although the Friedel oscillations contribute little to 
(n k ) , they are responsible for the fluctuations of the 
fermionic occupation numbers. Indeed, using Eqs. ([5]) 
and ITUl). we find 



(5n k Sn k ,) Q 



n 2 [(k - k')R] 
7r 2 (fc-fc') 2 



(13) 



Eq. (Tl3|) is valid when both \k\ and \k'\ are small com- 
pared with the Fermi momentum k F - In this limit the 
dependence of (5n k 6n k i) on k + k' [which we have ne- 
glected in writing Eq. (1 13[) ] is very weak. 
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The visibility of the Friedel oscillations Eq. (| 13|) is not 
affected by temperature T as long as T <C £f, where 
ep = Nlu — kp/2m is the Fermi energy (i.e., the chemi- 
cal potential for hard-core bosons). Indeed, a finite tem- 
perature introduces an uncertainty 5R ~ RT/ep in the 
size of the cloud which leads to the exponential decay 
of (8n k 5n k ) at \k - 



| > I /SR. The oscillations Ijl3|) 
survive as long as 8R <Si R, i.e., at T -C ep. 

So far, we demonstrated that the information about 
the broken translational invariance in the initial state of 
the system is preserved in the statistics of the fermionic 
occupation numbers. In particular, it manifests itself in 
the characteristic oscillatory dependence of (5n k 6n k i) on 
k — k', see Eq. (fl3|) . However, the fermions emerged 
in our problem merely as a convenient way of dealing 
with the exact eigenstates of the system. Since the rela- 
tion between the effective fermions and the original hard- 
core bosons is nonlocal, see Eq. ([3]), the behavior of the 
bosonic correlation functions is much more complex. 

We discuss here the bosonic occupation numbers f k = 
<^>\<P k - Unlike their fermionic counterparts Eqs. (JB]) and 
([H]). the expectation value (fk)t (momentum distribution) 
and the fluctuations (5f k 6f k >) t (here Sf k = fk - (fk)t) 
are no longer constant. However, the time-dependent 
contributions are superpositions of an infinite number of 
oscillating terms and decay at t —* oo [l2j • 

In order to find the bosonic occupations at t — > oo, we 
again concentrate on the harmonic trap potential Eq. ([9]). 
Following the method of 13, 14, lH, we consider first the 
single-particle Schrodinger equation 

d , , , 1 d 2 , , , 

with the initial condition ip n (x,0) = cf> n (x), where (j> n (x) 
is the normalized stationary eigenfunction of a harmonic 
oscillator corresponding to the eigenenergy e n = (n + 
l/2)u>. At t ^> 1/w, the wave function ip n (x,t) assumes 
the form 



ip„(x,t) 



1 



exp 



2l 2 uot 



e"' v i(i/wi) , (14) 



where r is given by Eq. (| 12[> . Upon introducing dimen- 
sionless variables 

x 

rj = ut, £ = — , 

m 

we rewrite Eq. (fl4|) as 

^ n (x, t) = rj-^e"'^ 2 e-^^l) . (15) 



Using Eq. (|15[) . the first-quantized many-particle wave 
function of hard-core bosons 17] &t can be expressed via 
its initial value 5 



(16) 



here {xi} — x\, . . . ,xn, £i = Xi/r]l, and &o({xi}) is the 
many-body eigenstate of the initial Hamiltonian H with 
energy E . 

In the second-quantized language, Eq. p6[) implies the 
operator relation 11 1 



(17) 



where (p{x, t) and (f>(x, t) are operators in the Heisenberg 
representation with the time dependence governed by the 
Hamiltonians Hq and H = Hq + V, respectively: 



(f{x,t) 



JH t 



<p(x) e 



-iHot 



iHl 



<p(x) e 



Mi 



Substitution of Eq. ([TTJ) into 



fk(t) - ^ j dxdx' e ik ^-^^{x,t)^{x' ,t) 



yields 



/*(*) 



2tt J 



(18) 



At r\ — > oo the integral over ^ and £' here can be evaluated 
in the stationary phase approximation with the result (TH 



f k (t^oo) = e^[l 2 p(kl 2 )]e-^ 



(19) 



[Analogous calculation for fermions yields Eq. (TT2]) which, 
unlike Eq. (fT9)l . is valid at all t > 0.] Finally, comparing 
Eq. ([Ill) with Eq. (p}, we find [ll| 



/*(* 



oo 



n k . 



(20) 



According to Eq. (f^Uj) . at t — > oo the bosonic occu- 
pation numbers in fc-space f k coincide with the integrals 
of motion n k . Since Eq. ([20]) holds for operators [ll| , it 
also implies that the statistical moments of the bosonic 
occupation numbers f k at t — > oo coincide with those for 
fermions in the initial trapped state, e.g., 



(/fe)oo = ( n k) ; (SfkSfk'}^ = (Sn k 5n k ,) . 



(21) 



We pause now to discuss the conditions of applicability 
of Eqs. ([2U)l and (|2"Tj) . Since we used the stationary 
phase method, the relevant characteristic time scales can 
be obtained by equating the scale of variation with £ of 
the phase ?y£ 2 in Eq. (fT7|) with that of the field (f(t;l, r). 

The dependence of ip{x) [and, therefore, of <p(x, r)] on 
x is characterized by two length scales. The longer one 
is the size of the trapped system R. Neglecting all other 
scales, we find t\ ~ 1/ep; this corresponds to the time 
the particles in the trap move between the collisions (the 
mean free time). The shorter scale of variation of ip(x) is 
the distance between particles R/N, which leads to the 
time scale ti ~ 1/oj ~ Nti. This is the time it takes for 
a particle moving with the Fermi velocity vp = kpjm 



4 



to cross the trap, t% <~ R/vf, or, equivalently, for the 
expanding cloud to double its size. 

Shortly after release of the trap, at t\ <C t -C t 2 , the 
oscillating transient contributions to the bosonic momen- 
tum distribution are still present, but (fk)t averaged over 
time is already given by the smooth "fcrmionic" semicir- 
cle Eq. (|lip . see 13j. In this regime the discreteness of 
the system is not important and the "shot noise" fluctu- 
ations, Eq. (fT3|) . are not yet resolved. Accordingly, the 
hydrodynamic description 13, 3] based on Eq. pro- 
vides a complete information about the system. 

Much later, at t 3> t 2 , the system enters the asymp- 
totic regime where Eq. (|20p is applicable. In this regime 
the transients have already decayed [T3| , the statistics of 
particles no longer matters, and statistical moments of 
the bosonic occupation numbers approach their steady- 
state values Eq. (f2~Tj) . In other words, this regime is es- 
sentially that of the collisionless expansion of the system. 

It should be noted that the setup discussed here is es- 
sentially the same as that studied recently in [r| • Based 
on the behavior of (/fc) t , it was conjectured there that 
any isolated system with integrable dynamics relaxes to 
a state described by a certain generalized Gibbs distri- 
bution. According to the prescription of [HI adopted 
for continuously varying k [19|, the density matrix at 
t — > oo has the form pa cx exp (— J dk /Jfenfc). Although 
the coefficients (3k can always be chosen in such a way 
that Tr(pGrifc) = (rifc) [18], finding (3k is obviously not 
sufficient to establish the validity of the conjecture. In- 
deed, no matter what (3k is, (flunk') = Tr(/3G^fe^fe') = 
(nfc)(nfc/), i.e., (5nk5nk) = 0. In view of Eqs. (Tj~3")) . 
(f^Uj) . and (jnj), we conclude that pq not only neglects 
the correlations between different integrals of motion rife, 
but also fails to account for the correlations between 
the bosonic occupation numbers fk- This raises serious 
doubts whether the generalized Gibbs distribution con- 
jectured in [l8| is actually useful for the description of 
the quantum quench problems. 

Although our consideration relied rather heavily on the 
properties of the hard-core Bose gas with the harmonic 
initial confinement, we expect some of our conclusions to 
be generic. In particular, we expect that any finite ID 
system with short-range interactions enters the collision- 
less expansion regime at t t 2 ~ R/v s , where v s is the 
sound velocity in the initial trapped state. (In an infinite 
system, this time scale corresponds to the establishment 
of a local equilibrium in a subsystem of size R Q). 

In a nonintegrable system, a relaxation would occur at 
t < t 2 . The relaxation would partially restore the trans- 
lational invariance, leading to the suppression of the noise 
{fifkdfk')oo- Since the noise is not sensitive to tempera- 
ture (see above), the accuracy of noise measurements in 
time-of-fiight experiments is not limited by one's ability 
to control the temperature. This suggests that deviations 
from the integrability are easier to detect in noise mea- 
surements than, for example, by observing the saturation 



of the height of the peak in the dynamic structure factor 
with loweringthe temperature [201 ] in Bragg spectroscopy 
experiments [4|. 

Finally, we emphasize that real-life ID bosons are nei- 
ther hard core nor their dynamics is integrable. It is 
conceivable that deviations from integrability will have a 
dramatic effect on the behavior of some observable quan- 
tities. Detailed understanding of the relaxation mecha- 
nisms and other consequences of nonintegrability in ID 
Bose systems remains an important open problem. 
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